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Abstract 



This paper is dedicated to the memory of our collaborator Nathan Isgur. 

Quark-hadron duality is an interesting and potentially very useful phe- 
nomenon, as it relates the properly averaged hadronic data to a perturbative 
QCD result in some kinematic regimes. While duality is well established ex- 
perimentally, our current theoretical understanding is still incomplete. We 
employ a simple model to qualitatively reproduce all the features of Bloom- 
Gilman duality as seen in electron scattering. In particular, we address the 
role of relativity, give an explicit analytic proof of the equality of the hadronic 
and partonic scaling curves, and show how the transition from coherent to in- 
coherent scattering takes place. 
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I. INTRODUCTION 



Quark-hadron duality has been well established experimentally for over 30 years, but 
our theoretical understanding of the phenomenon is quite limited so far. In the past year, 
there has been renewed interest in duality, both on the experimental and theoretical 
side Duality is a major point in the planned 12 GeV upgrade of CEBAF at Jefferson 

Lab 0. Duality between partons and hadrons is also employed in QCD sum rules. |10 . 

In a recent publication , we presented results obtained in a confining, relativistic model 
which qualitatively reproduced the features seen in Bloom-Gilman duality. In we only 
discussed a reaction where all particles involved - electrons, photons, and quarks - were 
treated as scalars. In this paper, we present the model for physical electrons and photons, 
and only treat the quarks as scalars. We describe the model and its properties in more detail, 
and discuss the Coulomb sum rule, the transition from coherent to incoherent scattering, 
and duality in the form factors. We also put special emphasis on the role of relativity. 
Relativistic treatment was one of four basic conditions imposed in to obtain duality, and 
here we show the consequences of relaxing this condition, and compare non-relativistic and 
relativistic calculations. One main point of this paper is the explicit derivation of the scaling 
curve for scattering from quarks confined in their initial and final state, and the comparison 
of this scaling curve to that obtained in a parton model calculation. 

For the convenience of the reader, we define the concept of duality in the following, 
and briefly discuss a few basic implications. In the literature, there exist many slightly 
varying "definitions" and usages of the term duality, and the phenomenon manifests itself 
experimentally in many different processes. We begin with a definition that covers all these 
cases. First, we need to make an obvious observation: any hadronic process can be correctly 
described in terms of quarks and gluons. In other words, quantum chromodynamics (QCD) 
is the correct theory for strong interactions. While this statement is obvious, it has little 
practical value, as in most cases, we cannot perform a full QCD calculation. E.g., in order 
to calculate a resonance excitation form factor, one would need to include very many quarks 
and gluons, and they would all couple strongly. We will refer to the above statement that 
any hadronic process can be described by a full QCD calculation as "degrees of freedom 
duality" . 

There exists a more practical and less obvious version of the first statement: in certain 
kinematic regions, the average of hadronic observables is described by a perturbative QCD 
(pQCD) result. This is the statement of duality, and we are going to explain the details in 
the following. 

With pQCD result, we indicate the result for the underlying quark process - for inclusive 
inelastic electron scattering from a proton, it is free electron-quark scattering; for semilep- 
tonic decays, e.g. B XcWi, it is the underlying quark decay rate, in this case obtained 
from the process b cWi |TI|]; for e~^e~ hadrons, it is the underlying e'^e~ — > qq process. 
Now, it is clear that we expect perturbative QCD to describe Nature in a certain kinematic 
regime, i.e. for very large Q^. In this regime, due to the fact that full QCD is approximated 
by perturbative QCD, the statement of duality turns into the statement of the "degrees of 
freedom duality" . So we have identified one kinematic regime in which even the non-obvious 
version of duality must hold. 

We also can identify a kinematic regime for which duality cannot hold: for 0. 
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While the underlying reason for the breakdown of duality at low four-momentum transfers 
is the non-perturbative, strong interaction of the hadrons, one can see the breakdown of 
duality most easily by considering the transition from incoherent to coherent scattering. For 
duality to hold for the nucleon structure functions in this case, we would need the following: 
the elastic proton and neutron form factors, which take the value of the nucleon charge for 

— > 0, would have to be reproduced by electron scattering off the corresponding u and d 
quarks. Now, for the proton this can work, as the squares of the charges of two u quarks and 
one d quark add up to 1. However, for the neutron, the squared quark charges cannot add 
up to 0, so it is clear that duality in inclusive inelastic electron scattering from a neutron 
must fail for 0. In addition, we know from gauge invariance that for 0, at fixed 

energy transfer u, the function vW2{i'-,Q'^) must approach 0. It is clear that the scaling 
function F2{x) does not show that behavior, which gives us an additional reason to expect 
the breakdown of duality at low Q^. 

So now we know that duality has to hold in one kinematic regime and that it has to 
break down in another kinematic regime. Obviously, a very interesting question is what 
happens in between these regimes, i.e. how exactly does duality break down, how far does it 
hold in the regime where it is nontrivial, i.e. for moderate values of Q^, and how accurately 
does it hold where it holds. 

II. THE MODEL 

Here, we present a model for the study of quark-hadron duality that uses only a few 
basic assumptions. Namely, we assume that it is necessary to include confinement and 
relativity in our model, that it is sufficient to base our model solely on valence quarks, and 
that these quarks can be treated as scalars. A model with these features will not give a 
realistic description of any data, but it should allow us to obtain duality and study the 
critical questions of when and how accurately duality holds. 

Although it is our aim to study duality in electron scattering from the nucleon, i.e. from 
a three-quark-system, as a first step we simplify the problem at hand by substituting two 
quarks by an antiquark, as the representation 3 ® 3 in SU(3) contains the representation 
3. This means we have a two-body problem now, and we have to solve the Bethe-Salpeter 
equation. In the special case of the mass of the antiquark, M, going to infinity, the problem 
further simplifies to a one-body problem. In the case of scalar quarks considered here, we 
obtain a Klein-Gordon equation. In contrast to where we assumed that all particles 
involved - electrons, photons and quarks - are scalars, here we treat only the valence quarks 
as scalars. Note that the experiment which our model resembles most would be electron 
scattering from a B meson. Still, we expect to gain valuable insights from considering this 
case, and would like to stress that none of the assumptions we made here prevent us from 
extending our model to describe more realistic circumstances. 

We have chosen to implement confinement by using a linear potential, which leads to a 
relativistic harmonic oscillator solution. This has the advantage that analytic solutions can 
be readily obtained and that a comparison with the non-relativistic case is easily feasible. 

We have to solve the Klein-Gordon equation with a scalar potential: 
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- V' + + V^^ ^x) = (1) 

with the usual ansatz ^^{x) — $(r ) exp(^iEt) and the confining potential 

V\f)=b\^, (2) 

where b is the relativistic string tension and has dimension [b] = [mass^]. The superscript 
of the wave function denotes positive and negative energy solutions to the Klein-Gordon 
equation. The mass of the quark is denoted by m and we use m = 0.33 GeV throughout 
this paper. 

The Klein-Gordon equation in this form can be easily rearranged to have the form of a 
Schrodinger equation to give 

V' lb' 2^ ^.^x E'-m^ 

2m 2m J ^ ^ 2m ^ " ^ ' 

where the similarity to the Schrodinger equation for a non-relativistic harmonic oscillator 
potential becomes apparent. The solutions to this equation are easily obtained by making 
the substitutions k = — and E = ^ and using the well known solutions of the non- 
relativistic case. 

The energy eigenvalues for the Klein-Gordon equation are E — ±En where 



EN^^2(3'{N+^)+m'. (4) 

N is the principal oscillator quantum number and P = b^/'^ . The corresponding wave func- 
tions are the usual nonrelativistic oscillator wave functions. For the present application, it 
is convenient to express the oscillator wave function in Cartesian form as 

*iv(r ) = 0„Jx) (t)n,{y) (pnM (5) 

where 

<P^^{x) = , ^ Hr^^iPx) exp {-l(3'x') (6) 
^j2^-njy/^ ^ 

with similar expressions for y and z coordinates. The iJ„ are the Hermitc polynomials. 
Unless noted otherwise, we use /3 = 0.4 GeV, which was chosen to give reasonable values for 
the mass splitting and charge radius. 

It should be noted that the negative energy solutions are just that since we are using 
a one-body wave equation and not a field theory. Therefore, these are an artifact of the 
model, but arc necessary to provide a complete set of relativistic states. 

As we choose to retain the non-relativistic wave functions, we differ from the usual 
relativistic normalization: using the Klein-Gordon normalization condition for these wave 
functions gives 
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This leads to the factor I/AEqEn in the response functions, and to the energy factor in 
the current operator defined below. Of course, we could have used explicitly relativistic 
normalized wave functions, which would have led to a different expression for the form 
factor, given below in Eq. (p!3[) . 

The basic difference between the relativistic and the non-relativistic oscillator equation is 
the difference in the energy spectrum: while the non-relativistic solutions are equally spaced, 
as oc iV, the relativistic spectrum goes as oc ^/N for large N so the density of 

states increases with increasing N. We note in passing that with this relativistic spectrum 
our simple model gives rise to linear Regge trajectories [T^ as seen in Nature. 



In the following, we will consider electron scattering from a meson with an infinitely 
heavy antiquark. In contrast to our previous publication where we treated all particles 
as scalars, the electrons in this paper are spin 1/2 fermions, and the virtual photons have 
spin 1. Unless otherwise noted, we assume in this paper that only the light quark carries 
a charge, and that the photon therefore couples only to the light quark, not to the heavy 
antiquark. 

For a photon coupling to the quark in the positive energy harmonic oscillator ground 
state and transferring the four-momentum q = (z/, g) , we have the following current matrix 
element: 

nq) = I d'x expi-tq ■ x) (<l>^/a;) 9^$o(x) - (S'^^^/x)) ^ix] 

where Nf can designate either a positive or negative energy state. Using this definition of 
the current along with (|l|), it can be easily shown that the current is conserved, g^j'^ = 0. 

The calculation of the double differential cross section is straightforward and leads to 
the Rosenbluth equation 

o-Mott [vLRL{q, jy) + vtRriq, j^)] , (9) 



dEfdQf 

where (JMott is the Mott cross section, vl and vt are the usual leptonic coefficients 

q^ 2 q^ 2 

= —q^ = g*^ — z/^, and the longitudinal and transverse response functions are 

oo -| 

Rdq. ^)=Y. , r, r, l^oNm' liEo + E^f5{p + Eo- E^) - (Eo - E^f5{p + Eo + E 



(10) 



and 



a °° 1 

Rriq, = 8^ E T-^r-ET- N |Foiv(gir [S{iy + - E^) - 6{u + Eo + E^)] . (11) 
q N=o 4 Eq En 

In these expressions, Fq^n stands for the excitation form factor, 

FoAq') = I d'fexpizq- r) ^%^{r) %{f) = J d^p^%^,{p) ^o{p- q) , (12) 
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where \I/ indicates a momentum space wave function. Making use of the recurrence relations 
of the Hermite polynomials, we find an explicit expression for Fq^j^: 



F„,.fe-)^^."(J^)"exp(-il). (13) 

Note that some care is necessary in writing the expressions for the responses to properly 
include the negative energy states. The relative sign between the positive and negative 
energy contributions is associated with the negative norm of the negative energy states. 

These expressions for the response functions have been derived assuming that the quark 
is excited from the ground state into a resonance state, N, and remains there without 
decaying. This is just the first step on the way to meson production in this picture. The 
^-function in the energies is an artifact of this assumption. 

Note that because we assume scalar quarks, there is no magnetization current present. 
The only contribution to the transverse part of the cross section comes from the convection 
current. As a result, the transverse response falls faster than the longitudinal response with 
increasing momentum transfer, as will be shown explicitly below. This is in contrast to the 
case of spin-1/2 quarks where the magnetization current dominates. In turn, it causes the 
transverse response to dominate at large momentum transfer, giving rise to the Callan-Gross 



equation ||T3[ in the scaling region. 

The inclusive, inelastic electron scattering cross section can be re-expressed in terms of 
two structure functions, Wi and W2, which depend only on z/ and Q^: 

"^"^ cTMott (W2iu, + 2 W.iu, tan^ ^ ) . (14) 



where 



W,{u,Q') = ^RTi^jQ^ + u^,u) 



III. THE COULOMB SUM RULE 

For the moment, we will consider a wider class of models for hadrons made up from 
confined quarks, namely the more general case of models where all quarks carry an electric 
charge. It is interesting to consider an apparent contradiction between a model such as the 
one discussed here, and the parton model. In our model, since all states are bound states, all 
of the transition form factors are coherent in that they are the result of scattering from the 
total charge. The parton model however assumes that the cross sections are composed of 
incoherent scattering from the individual constituents resulting in cross sections proportional 
to the sum of squares of individual charges. One method of examining the transition from 



coherent to incoherent scattering is the Coulomb Sum Rule |TJ]. Consider the longitudinal 
response function 
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RL{q, iy)^J2< ^o\pKQ)\'^f >< i^f\p{Q)\^o > + Eo - Ef) (16) 
/ 

where the sum represents a generahzed sum over all final states, bound or free, and p{q) is 
the Fourier transform of the charge operator. Now define the logitudinal sum as 

/oo 
duRL{,q,p) (17) 
-oo 

Using the above definition of the longitudinal response and the completeness of the final 
states, this becomes 

/oo 
dvY,< ^o\p\mf X i^flPmi^o > + Eo- Ef) 
f 

= <V'o|p^(g>(?)|V'o> ■ (18) 

This is a general result. To see how this relates to the problem of the transition between 
coherent and incoherent scattering, consider the simple case of a nonrelativistic system of 
two constituents with charges ei and 62. In this case 

S{q) = J dh^l{r) (eie-'^*-''! + e2e-''l-'^) (e^e^'l-''' + e2e"^-''^) M^) 

= I d\4{r) [el + el + eie2e-'^<'''-'''^ + eie2e'^-^''^-'''^) Mr) 

= el + el + 2eie2 J^{q) , (19) 

where f — fi — f2 is the relative coordinate of the two particles and 

J^{q) = ^J d^rijlif) e'^-^M^ (20) 

is the real part of the Fourier transform of the ground-state probability density. 

In order to understand the physical significance of the quantity determining the rate of 
fall-off of the mixed term containing the product of ei and 62, it is necessary to write down 
the most general form of the charge form factor for two quarks with charges ei, 62 and masses 
mi,m2 (M — mi + 777.2). Here, we have dropped the (5-function obtained from integrating 
over the cm. motion in the second step: 

Fo,nM = / d'f'^h^'^ (eie^«'-^^i+e2e^^"^^^)*o(r) 

^ J d¥*^/r) (eie^^"^ + 626"^^"^) *o(r) (21) 

From this expression one sees that in the most general case, J-'{q) cannot be interpreted in 
terms of the ground-state charge form factor. However, in the special case of mi = m2, 

FoAq) = (ei + 62) / dVe^^--5>o(OI' = i^i + e2)hM ■ (22) 
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In this case 



■151 



Hq) = /o,o(2g1 . (23) 

A special case of our general result ( [T^ ) was discussed in where the case of two scalar, 
equal mass quarks in a non-relativistic harmonic oscillator potential was considered. 

In the model we present in this paper, all of the charge is carried by one of the con- 
stituents, that is ei = 1 and 62 = 0. Since there is only a single charge, there is no difference 
between coherent and incoherent scattering, so we expect that 

S{q} = 1 . (24) 



This then provides a useful test of the model. 
Using (|10D, 

/oo °° 1 



X [{Eo + ENf5{v + Eo- En) - {Eq - E^fSiu + Eq + E^) 

^ {Eq + -Ejv)^ - (Eq - EnY I ^ ,^.2 
1^ AW p FOTVlgJl 



00 



= E \FoN{q)? (25) 

N=0 

Using ( p!3D it is straightforward to demonstrate that S{q) = 1 in this case so the Coulomb 
Sum Rule is satisfied. Indeed, this will be true regardless of the form of the confining 
potential, as long as one considers a complete set of solutions. Note that for this model it 
is necessary that the integral in ([17|) be over both positive and negative energy transfers for 
the sum rule to be satisfied. 

In electron scattering, only the spacelike region is accessible and the negative energy 
states are an artifact of the use of the Klein-Gordon equation as a wave equation. It is 
useful, therefore, to examine the contributions to the Coulomb sum from spacelike, timelike 
and negative energy states. When referring to the spacelike and timelike contributions, 
only positive energy states are included. The different contributions are shown in Fig. 
At |g| = 0, only the elastic form factor can contribute and therefore must saturate the 
sum rule. As the momentum transfer increases the elastic form factor falls off, resulting in a 
decrease in the spacelike contribution. It decreases until the momentum transfer increases to 
a point where the first excited state enters the spacelike region. The spacelike contribution 
then saturates the sum rule again. This process continues as new form factors become 
accessible in the spacelike region. The result is a saw-toothed behavior of the spacelike 
contribution. Because the density of states for this oscillator model increases with increasing 
energy, the magnitude of the "teeth" becomes smaller with increasing momentum transfer 
until the spacelike contribution is essentially smooth. Note that the spacelike contribution 
over-saturates the sum rule around momentum transfers of 0.5 GeV and above. Since the 
jaggedness in the spacelike region is associated with the migration of contributions from 
the timelike region, it is not surprising to see that the complement of this behavior does 
indeed show up in the timelike contribution. As the momentum transfer increases, the size 
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FIG. 1. Positive energy spacelike, timelike and negative energy contributions to the Coulomb 
Sum rule from our model. 



of the timelike contribution becomes small. The negative energy contribution is smooth and 
compensates for the over-saturation of the spacelike contribution. This is clearly an artifact 
of using a one-body wave equation with negative energy contributions. 

It should be pointed out here that the small size of the timelike contributions is an es- 
sential consequence of using a relativistic wave equation. This can be seen by considering 
a similar situation where the nonrelativistic oscillator model is used to describe the sys- 
tem. The spacelike and timelike contributions of such a model are shown in Fig. ^ The 
Schrodinger equation of course has only positive energy solutions. Here the saw-toothed 
behavior in the spacelike contributions is similar to the relativistic case with the important 
difference that due to the linear character of the nonrelativistic spectrum (see Fig. the 
contributions from states entering the spacelike region is not rapid enough to compensate 
for the fall-off in the form factors. Therefore, with increasing momentum transfer, the size of 
the spacelike contribution approaches zero with all of the strength appearing in the timelike 
region. 

Duality in the sum of the form factors in the spacelike region is related to the Coulomb 
sum rule. While the concept of duality in the form factors is not clearly related to an 
observable (in contrast to duality in the structure functions), it still has received attention 
in the literature ||^. We describe it within our model in appendix]^. 

This model can be easily extended to the case where both quarks are charged by exam- 
ining the behavior of the two-body Gross equation in the limit where the mass of one of the 
particles becomes infinite [|16|. The contribution of the infinite mass particle to the struc- 



ture function is simple and straight forward. Due to its infinite mass, this particle remains 
stationary for any finite momentum transfer, and it is point-like. This particle therefore 
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FIG. 2. Spacelike and timelike contributions to the Coulomb Sum Rule for a nonrelativistic 
oscillator. 



contributes only to elastic scattering and has a constant form factor. The structure function 
can then be written as 



Riiq^v) = \eiFoo{q) +62^ 
00 1 

N=l 



- (Eo - E^fSiu + E0 + E 



N, 



The Coulomb Sum can be easily calculated to be 



S{q) =elY: \FoN{q}\' + el + 2e^e2^FoM 

N=0 

= el + el + 2eie2^Foo{q) 



So in this case 



(26) 



(27) 



^(g) = moiq) = Fooiq) . 



(28) 



After examining how the apparent contradiction between coherent and incoherent scattering 
is resolved in a more general framework, we now proceed to investigate duality in our model. 
The first condition for duality is that one obtains scaling in the structure function calculated 
solely with resonances, and that the scaling curve thus obtained agrees with the scaling curve 
obtained in the parton model. 
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IV. THE PARTON MODEL 



The usual assumption of the parton model is that at large momentum transfers the final 
state quarks can be treated as though they were free. Examination of the structure functions 
for large and fixed Bjorken x leads to identification of the scaling functions. 

For our simple model, the response functions for excitation of a bound (off-mass-shell) 
quark to a plane-wave final state can be calculated analytically as 



exp 



X [e{v + Eo-m)- 9{-u -Eq- m)\ 



r 



(29) 



and 



/5 



2n-2Eo\q\^ [ 



- \ 2{y + \q\) + '—j exp 



/ y2 



X [e{iy + Eo-m)- e{-v -Eo- m)] 



(30) 



where y — ^Jii^ -\- EqY — m? — \q\. 

For our model it is not possible to define the scaling variable in terms of the target mass 
since it is infinite. For this reason we define a new Bjorken variable 



M 



Q2 

UBj = XBj = 

m zmu 



which covers the interval — oo < < oo. Using v — , and taking the limit 
the structure functions become 



and 

RT{Q^,UBj) 

Since in this limit 



87T2 (3mEoUBj 
2(5muBj 



exp 



{Eq - muBjf 

/?2 



'2mu 



Bj 



2mu 



Bj 



(31) 



oo, 



(32) 



exp 



{Eq - muBjY" 



n-^EoQ^ 



P' 



0(7^) -0{- 



'2mu 



and — 



Bj 



m\UBj\ 



2mu 



Bj J 



(33) 



the structure functions have the limits 

Wi{uBj,Q^)^ F^{uBj)^Q 

and 



(34) 



(35) 
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u\W2{u,Q^)^F2{uBj) 




{Eq - muBjf 

/32 



(36) 



Note that the choice of \v\ in defining W2 is necessary to provide a properly normahzed 
scahng function as will be seen below. 

Although we have used the Bjorken scaling variable to obtain these results, this will be 
true for all such variables since all acceptable scaling variables must reduce to the Bjorken 
scaling variable as — >^ 00. Therefore, a more general expression for F2 for any scaling 
variable and any initial state can be written as 



After obtaining the scaling curve in the parton model, i.e. the scaling curve for a quark 
initially bound and then free, we proceed to find an expression for the scaling curve in our 
model, where the quark makes the transition from the ground stated to an excited bound 
state. 



An interesting feature of our relativistic oscillator model is that the scaling behavior of 
the model can be determined analytically by making a continuum approximation. The jus- 
tification for this is that at increasing momentum transfer the contributions to the response 
functions are dominated by higher energy states. Since the density of states increases with 
increasing energy, it is reasonable that a continuum approximation should provide a good 
description of the averaged response for large momentum transfers. 

Using (|10D and ([T^) we can write 




(37) 



where N{p) is the ground state momentum distribution normalized such that 




(38) 



V. CONTINUUM LIMIT 



00 ^ 1 / g 2 \ -2 



X [{En + Eof5{EN - Eo - u) - {Eo - EnYS^En + ^0 + i^) 



(39) 



where AiV = 1 



It is convenient to write 




(40) 



where 



(41) 



^From this it can be determined that for a variation in A^, 
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and 



AE 



N 



Pn^Pn 



E 



N 



The longitudinal response function can then be rewritten as 



E 



AE„ 



-2 \ 



exp(-^) 



X [(En + EofSiEN - E^ - v) - (Eo - En)''S(En + Eo + u) 
This sum can now be approximated by the integral 



AP'^Ei 



roo 

■0 Jeo r(i 



dE 



(l + ^ 



2/32 



X 



[E + E^fbiE -Eq-v)-{E- E^fbiE + Eo + I/) 
which can be trivially evaluated to give 



{v^'lE^f{Q{v) 



-2Eo - v)] ( q 



(u+Eq)'^-eI 
■2 \ 5^2 



q- 



4/32£;or (i + W 

Similarly, the transverse response becomes 



exp(-^). 



+ 2Ef,v) [e{v) - e{-2Eo - 1^)] ( r 



20^ 



Eoq-^r (l + (^^±^) V2/5^ 



exp(-^) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



In the scaling limit the argument of the F function becomes large, so Stirling's formula 
can be used to write the longitudinal response function as 



IEq) «Xp ^^ 2/32 "-^ \v-2+2Eov ) 



2/3- 



2/32 



X \Q{v) - ei-2Eo - y)] . 



Using 



u — 



2muBj ' 



(48) 



(49) 



and taking the limit — > 00, the structure functions become 
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Similarly, 



8tc2 pmEoUBj 



exp 



2l3muBj 



exp 



(£"0 - muBj 



{Eq - muBj] 
(3^ 



2mu 



Bj 



2mu 



Bj 



(50) 



0{ 



2mu 



Bj 



2mu 



Bj J 



(51) 



Since (|50[) and ( pi]) are identical to (^) and (|33D the model scales to the parton model result. 
So, even though our model describes a bound quark being excited to resonance states, we do 
obtain a scaling curve in the Bjorken limit. This, as well as the results presented in |T7|-p!9|,|6| , 
show that scaling does not necessarily imply scattering off free constituents, a belief which 
is encountered widely. 

While others [|T7|-[T9[| have studied the transition from ground state to excited bound 
states and found scaling in similar models. However, duality is fulfilled only when i) the 
transition from ground state to excited bound states scales, ii) the transition from ground 
state to a plane-wave final state scales, and iii) both scaling curves coincide. In this section 
and the preceding section, we have shown that duality holds explicitly in our model. The 
numerical approach towards the scaling curve is shown in Fig. |^. 

Note that as we explicitly made use of Stirling's formula in the derivation of the scaling 
function in the continuum limit, it is clear that for the lower lying resonances, which corre- 
spond to lower values, we will never quite see scaling in the subasymptotic regime. This 
is of no practical relevance, as these resonances are pushed out to very high values of u for 
larger Q^, and the structure functions practically vanish in this region. This is completely 
analogous to the fact that for electron scattering from a proton, one always picks up the 
elastic scattering at a;^^ = 1, independent of - even though the elastic form factor will 
have fallen off to negligible values at high enough Q^. 



VI. APPROACH TO SCALING IN THE STRUCTURE FUNCTIONS 

After establishing analytically that one of the necessary conditions for duality is fulfilled, 
namely scaling to the scaling curve obtained from a free quark in the final state, we proceed 
to investigate the approach to scaling numerically. 

We would like to remind the reader that our results should not be compared to the 
available nucleon data - our model calculations describe a process that might resemble 
electron scattering from a B meson, which has never been measured. In general, when we 
consider scattering from a meson target, scaling will set in later than for a baryon target: 
momentum sharing for higher momenta between fewer constituents is easier, which leads 
to a slower fall-off of the individual form factors, and to a later onset of scaling. In our 
situation, where the system is not allowed to decay, we have a somewhat extreme case. 

To see duality clearly both experimentally and theoretically, one needs to go beyond 
the Bjorken scaling variable xbj and the scaling function S2Bj = that goes with it. 
This is because in deriving Bjorken's variable and scaling function, one not only assumes 
to be larger than any mass scale in the problem, but also that high (pQCD) dynamics 
controls the interactions. However, duality has its onset in the region of low to moderate Q^, 
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and there masses and violations of asymptotic freedom do play a role. Bloom and Gilman 
used a new, ad hoc scaling variable u' ^ in an attempt to deal with this fact. In most 
contemporary data analyses, the Nachtmann variable [0,^ is used together with S2Bj- 



Nachtmann's variable contains the target mass as a scale, but neglects quark masses. For 
our model, the constituent quark mass (assumed to result of spontaneous chiral 

symmetry breaking) is vital at low energy, and a scaling variable that does not make any 
assumptions about the size of the quark and target masses compared to is desirable. Such 
a variable was derived more than twenty years ago by Barbieri et al. to take into account 



the masses of heavy quarks; we use it here given that after spontaneous chiral symmetry 
breaking the nearly massless light quarks have become massive constituent quarks, calling 

it Xr 



l-2i?(v/^^^W + ^)(l + v/rT^) forz.<0 • 

where the definition for negative energy is chosen such that it satisfies the kinematic con- 
straints in this region and reproduces the behavior of xbj for large Q^. The scaling function 
associated with this variable is given by: 



S2cq=\q\W2 = ,Ju^ + Q^W2. (53) 

This scaling function and variable were derived for scalar quarks which are free, but have 
a momentum distribution. The derivation of a new scaling variable and function for bound 
quarks will be published elsewhere. Numerically, this scaling variable does not differ very 
much from the one in Eq. (^2|). Of course all versions of the scaling variable must converge 
to XBj and all versions of the scaling function must converge towards Sbj for large enough 
Q^. One can also easily verify that in the limit m — one obtains from ([5^ ) the Nachtmann 
scaling variable. In the following, we use the variable Xcq and the scaling function S2cq- 

We are now ready to look at scaling and duality in our model. Since the target has mass 
M oo, it is convenient to rescale the scaling variable Xcq by a factor M/m: 

M , . 

U = Xcq . (54) 

m 

The variable u takes values from to a maximal, dependent value, which can go to 
infinity. The high energy scaling behavior of the appropriately rescaled structure function 
Scq is illustrated in Fig. ||. The structure function has been evaluated using the phenomeno- 
logically reasonable parameters m = 0.33 GeV and /5 = 0.4 GeV. To display it in a visually 
meaningful manner, the energy-dependent 5-function has been smoothed out by introducing 
an unphysical Breit-Wigner shape with an arbitrary but small width, F, chosen for purposes 
of illustration: 



where the factor / = vr/[^ + arctan ^^^^^jr-^] ensures that the integral over the (5-function 
is identical to that over the Breit-Wigner shape. As for the all scalar case discussed in [Q, 
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FIG. 3. The high energy scaUng behavior of S2cq as a function of u for various values of Q^. In 
the left panel, we have used T = 100 MeV to give the impression of real resonances even though 
this large value distorts the scaling curve somewhat. For any width equal to or smaller than this, 
the distortion is rather innocuous, and for F — > 0, the structure function approaches the scaling 
function (solid line) in Eq. p^), as shown in the right hand panel. The structure functions in the 
left panel are shown for = 100 GeV^. 

the smearing out of the 5-function in energy with the Breit-Wigner shape leads to a slight 
widening of the curve and flattening of the peak height. However, when choosing a smaller 
value for the Breit-Wigner width, these effects disappear, as seen in the right panel of Fig. |^. 

Not unexpectedly, the scaling curve we find here, when using photons and electrons with 
their appropriate spins, differs from the one in the all scalar case both in its final shape 
and the approach to scaling. Now, we have two terms, the longitudinal and transverse 
response function, contributing to the structure function W2 and therefore to the scaling 
curve. More importantly, the terms themselves are different and more complicated in the 
case considered here. The longitudinal part of the structure function contains an additional 
factor [Eq + En^, which was not present in the all scalar case. As shown in section 0, the 
transverse response vanishes like 1/Q^ in the limit 00, while the longitudinal response 

grows like Q^. This leads to a vanishing of Wi, and to a independent value for |g|VI^2- 
Even though, from eqs. ([T5|j50|j5T|) it is clear that at lower Q^, the transverse contribution 
to W2 will not vanish immediately, therefore making the approach to scaling slower. The 
effect is rather significant, though, as the contributions of the convection current are very 
small. They delay scaling slightly for the low u part of the curve. The effect would be more 
important for contributions of similar size within a certain kinematic range. In our case, 
the main effect of the transverse contribution is to slightly broaden the curve. For smaller 
Q^, this effect is more pronounced for low values of the scaling variable u, as the higher u 
correspond to lower lying resonances, which have only tiny contributions from the transverse 
part. 

As already mentioned above, for a proton target, the dominant contribution to the 
transverse response and overall is the magnetization current, which does not contribute for 
our scalar "quarks" . Note that both the transverse and longitudinal contribution to W2 are 
positive definite. If a dominant contribution in the transverse response is present, it should 
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lead to a different scaling behavior in the structure function Wi than in z/W2. For i'W2, the 
longitudinal term with different behavior will most likely slow the approach to scaling 
down, as it is going to be of comparable size to the magnetization current contribution at 
low Q^. This is a completely general observation, and one would expect to see faster scaling 
in Fi once the data are available. The same conclusion was reached on a different basis in 

i- 

The shape of the scaling curve is also different than for the all scalar case. The peak is 
higher, the curve extends to larger values of the scaling variable, and for u — > 0, the scaling 
function actually vanishes now, as expected from a valence quark distribution, even though 
we do not find a behavior oc y/u as seen for proton targets. However, we cannot expect 
to reproduce the correct distribution function for quarks in our simple model with scalar 
"quarks" . 

From the explicit expression for the scaling curve, it is clear that it depends both on 
the binding strength j3 of the harmonic oscillator, and on the quark mass. It peaks at 
Upeak = (Eq + \J El + (3^) /2m, slightly above the value u^l^J;^^ = ^ which we found for the 
all scalar case. Naively, for a target of mass M made up of non-interacting quarks of mass 
niq, one expects a spike at xbj = jf- In our case, the role of the mass of the quark 
is played by the ground state energy Eq, which appears everywhere (e.g. in fiux factors, 
normalization) where one would have the mass in the free particle case. As our variable is 
rescaled with the factor ^, we expect Upeak ~ It is interesting to note that this value 
receives a slight binding correction due to the conserved current employed here. We note 
that for weaker binding, the peak gets narrower and its position slides towards u ^ 1. As 
expected, in the limit of a free particle, (3 —>■ 0, we do obtain a spike of infinite height at 
M = 1 [in this limit, the scaling function becomes 6{u)/{mEQ)]. 



VII. MOMENTS AND FURTHER SUM RULES 

Now, we will discuss global duality, where the term global implies that we consider an 
average/integral over many resonances, which is compared with the corresponding integral 
over the scaling curve. 

Local duality implies that we compare the contribution of one single resonance or just a 
few resonances with the scaling result, i.e. with the free quark result. This will be discussed 



in section [VIII| . The concept of local duality is taken to its extreme when one focuses not just 
on one single resonance, but on one point only of the contribution of the single resonance, as 
it was done first by Bloom and Gilman when they compared the peak value of a single 
resonance with the value of the scaling curve. This version of duality has been investigated 



in Ref. |2^. As this ratio would depend strongly on the Breit-Wigner width we use to 
smooth out the 5-functions, it is not appropriate to consider it in this paper. 

Global duality was first quantified by Bloom and Gilman in the form of finite energy 
sum rules, where the integral over the scaling curve was compared to the integral over the 
resonance contribution. The integration range in both cases comprises the region of the 
scaling variable u' or u, respectively, which corresponds to the resonance region, defined as 
having an invariant mass W < 2 GeV: 
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^' duyW2{v,Q^)= duj'uW2{uj'), (56) 



7o 

where Wm ~ 2 GeV, and Vm = *" 2Mj^ • Here, denotes the mass of the nucleon 

target. The agreement between the left and right hand side of this equation is better than 
10 %; for the larger values of Q^, starting around ^ 2 GeV^, the agreement is quite 
impressive: 2% or better. 

While it certainly would be desirable to calculate the same finite energy sum rule in our 
model, there is a practical problem and a philosophical problem. Firstly, in our model, we 
deal with an infinitely heavy system, so that in principle, the invariant mass W of the final 
state is always infinity. Even if we could define a reasonable substitute for the invariant mass 
of the final state, picking an integration limit is a problem in principle: for our model, the 
scaling curve consists solely of resonance contributions, even though they cannot be resolved 
and form a smooth curve. So, any distinction between "resonance region" and "continuum" 
or "multiparticle final states" is artificial. The conventional definition of "resonance region" 
as the region where W < 2 GeV means the region where the resonances are prominent 
and dominant. However, it does not mean that for W > 2 GeV, there are no resonances 
present, and it also does not mean that for W < 2 GeV, i.e. in the resonance region, there 
are no non-resonant contributions at all. Experimentally, there is background for W < 2 
GeV, and there are resonances for W > 2 GeV, e.g. several A^* and A resonances. Also 
from a theoretical point of view, it is obvious that resonances which decay by creation of 
a quark-antiquark pair in the final state must be accompanied by a corresponding non- 
resonant production mechanism, where the pair creation takes place in the initial state, and 
the photon interacts with the preformed meson. 

Since the distinction between a "resonance region" and a "continuum region" has its 
problems, we utilize the moments of the scaling function Scq, where the integration range 
comprises the whole interval of the scaling variable. The physical information contained in 
the finite energy sum rule and the moments M„ is the same. The moments of a scaling 
function S{Q^, x) with a scaling variable x are defined as 

Mr,{Q^)= dxx^-^S2{Q\x). (57) 



Here, in contrast to [^, we do not include the unphysical region ]xmax, 1] in the integration 
interval. It is obvious that higher moments, i.e. n = 4, 6, . . ., tend to emphasize the resonance 
region, as for fixed Q^, the resonances are found at large x. The values of the moments 
decrease with increasing n. In our case, we change to u-type scaling variables, see Eq. |5^, 
so that 

Mr"(g') = du S2{Q\ u) , (58) 

Jo 

where Mmax corresponds to the maximum value of u which is kinematically accessible at a 
given Q^. By changing from x io u scaling variables, we change the upper integration limit 
from a value equal to or lower than 1 to a value considerably larger than 1 for > 1 — 2 
GeV^. This means that our higher moments will emphasize the low- lying resonances even 
more than the conventional, x-based moments. Also, the higher moments will be larger than 
the moments with small n. 
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Evaluating the moments of the structure function (^) exphcitly one has 
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The elastic contribution is given by 
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(59) 



(60) 



Note that for vanishing four-momentum transfer Q^, all moments take the value 1, indepen- 
dent of n. 




FIG. 4. The lowest moments M„ function of Q . 

In Fig. H, we show the moments for n = 2,4,6,8, which were obtained by integrating 
over the positive energy states only. One can see clearly that all moments flatten out, even 
though they did not quite reach their asymptotic value at the highest value shown. The 
lowest moment M„, is within 9 % of its asymptotic value at = 5 GeV^, and within 
% of its asymptotic value at = 20 GeV^. As expected, the higher moments, which by 
construction get more contributions from the lower lying resonances, need higher values 
in order to reach their asymptotic values. For Me, we find that it has reached 64 % of its 
asymptotic value at = 5 GeV^, and 88 % of its asymptotic value at = 20 GeV^. 



19 



From these numbers, we can see that even though scahng does not set in for < 50 GeV^, 
the asymptotic values at least of the lower moments are reached much earlier. This reflects 
the fact that Scq{u, Q^) approaches the scaling curve by shifting towards higher m, not by 
approaching it from below or above. 

Since the continuum approximation provides a relatively simple analytic expression for 
the structure functions, it is possible to use this to study certain properties of the moments. 
First, however, it is necessary to determine the validity of this approximation for the calcula- 
tion of moments. Figure ^ shows calculations of the first three moments Mq, Mi and M2. In 
each panel the solid line represents the exact calculation according to (^). The dash-dotted 
curve is a calculation of the continuum approximation with both longitudinal and transverse 
contributions, while the dashed curve includes only the longitudinal contribution. Note that 
the continuum approximation works very well down to a couple of GeV^. Note also that 
while the inclusion of the transverse contribution slows convergence to the asymptotic value 
for Mq it improves convergence for the higher moments. 

The continuum approximation can then be used to obtain an expansion of the moments 
in powers of l/Q^ reminiscent of the operator product expansion (OPE) series, 

Mi'm = j:%- (61) 

Note that we do not have any gluons in our model, and thus no radiative corrections. The 
expansion coefficients correspond to the non-perturbative matrix elements of higher twist 
operators in the OPE. Since this is an asymptotic series, the expansion will fail at low 
with the point at which the series diverges being dependent upon the order of the series. 
The expansion coefficients for the five lowest moments are shown in Table |. Contributions 
to the coefficients from transverse and longitudinal responses are shown along with the total. 
The obvious feature of these coefficients is that they are not in general small nor do they 
show any obvious convergence. The reason for this can be seen from Fig. ^. Here the exact 
result is shown as a solid line and is compared to the expansion with from one to four terms 
for both Mq and M4. Since the moments all have a finite value at = 0, the function 
cannot be analytic in Any expansion in this variable to a finite number of terms will 

at some point diverge, either above or below the correct result. Using an additional term to 
extend the approximation to lower must require that the coefficient of this term be of 
opposite sign to the preceding term, leading to an alternating series. Higher moments have 
the curvature toward the finite result occurring at increasing values of Q^. This requires that 
the size of the coefficients for the higher terms in the series must also be increasing. This 
shows that the global duality observed in our model is the result of a delicate cancellation 
between many "higher twist" terms. However, while it is fascinating to speculate if duality 
in Nature is realized by small higher order expansion coefficients or by cancellations, our 
model is too simple to allow us to draw any conclusions about this. In fact, results presented 
in p5[ indicate that the expansion coefficients have the same sign. One may hope that the 
building of more realistic models will allow us to gain a better insight into this question in 
the future. 

The relation asymptotic behavior of the moments to sum rules can also be addressed in 
this model. Consider the moments of F2{u) 
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FIG. 5. The first three moments, Mq 
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Ml, and M2. The sohd hne shows the exact result cal- 
culated according to Eq. 59. The dash-dotted lines show the moment calculated in the continuum 
approximation, the dotted curves show the purely longitudinal contribution to the moment in the 
continuum approximation. Inclusion of the transverse contributions improves convergence of the 
moments to their asymptotic values. 21 
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where the integral starts at —00 to include contributions from negative energy states as in 
the Coulomb Sum Rule. Using ( |37D this becomes 
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are proportional to the normalization integral of the momentum distribution. Comparing 
these to the corresponding values of Cg and c\ in Table | shows that the contributions from 
negative energy states are small. Note also that for a spin-1/2 constituent the expression 
corresponding to (|37| ) have a leading factor of u rather than v? as in this case. So the sum 
rules would be associated with Mi and M2 as expected. 
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VIII. THE LOW Q2 region 



After studying the scaling behavior of our model at high and the moments over a 
range of four- momentum transfers, we now discuss the behavior at low Q^. In this region, 
resonances are dominant for a wide range in the scaling variable. 

Before discussing the numerical results, a remark on the kinematics is in order. For a 
fixed resonance in inclusive electron scattering from the nucleon, its position in terms of 
Bjorken's scaling variable is given by Xj-es = wa 102 • This means that for higher Q^, 
the resonance position moves towards higher values of XBj-, and for very large Q^, xbj —>■ 1- 
In our case, the maximal value of the scaling variable u is larger than 1, and for very large 
Q2, the resonances move out to very large values of m, where their contribution is extremely 
small. 

If local duality holds, we expect the resonance curve to oscillate around the scaling curve 
and to average to it, once is large enough. In Bloom-Gilman duality, the finite energy 
sum rule gets within 5 % for > 1.75 GeV^. For lower Q^, the resonances approach the 
scaling curve from below. In our case, we have the onset of scaling for larger values of 
than observed in Nature. This is not unexpected, as we consider an infinitely heavy meson 
as target, and assume that this meson is made up of scalar quarks. For this reason, our 
cross section for photon exchange for large is dominated by the longitudinal part, and 
the transverse part, comprising solely the convection current, is very small. In Nature, for 
spin 1/2 quarks, we have the magnetization current, which is the dominant component of 
the cross section, and which therefore determines the scaling behavior. So we cannot expect 
our model calculation to show the same behavior as experimental data for the same values 
of the four-momentum transfer. 

In discussing local duality and resonances, the smoothing method used becomes impor- 
tant. The visual appearance of "resonances" depends on the chosen smoothing method: 
a bumpy structure is seen only when a Breit-Wigner shape is inserted for the energy 5- 
function. It also depends on the width chosen in the Breit-Wigner smoothing method. For 
a smaller width, the resonances are visible for higher Q^. Depending on the width of the 
Breit-Wigner, e.g. for F = 100 MeV, we do not see any resonances for = 5 GeV^, even 
though this value of is below the scaling region. In this paper, the working definition 
of local duality which we use is "resonance curves oscillating around the scaling curve". 
At some point, when considering more realistic models, it may be useful and necessary to 
introduce a sharper, more quantitative definition. However, at this stage, we are interested 
more in qualitative results, and do not intend to quantify how well exactly local duality 
works for our simple model. 

In Fig. 0, we show our results for the scaling function S2cq{u, Q'^) for various low values 
of Q^. The (5-function in the energy has been smoothed out using the Breit-Wigner method, 
with a width of F = 100 MeV. For visual purposes, we have assigned a small width to the 
elastic peak, too. One can see clearly from the figure that the resonances move out towards 
higher u with increasing Q^, as dictated by kinematics. While the elastic peak is rather 
prominent for = 0.5 GeV^ and = 1.0 GeV^, it becomes neghgible for > 2.0 GeV^. 
This is the phenomenon we have observed already when studying the moments: the elastic 
contribution there vanishes rapidly with increasing Q^. 

As already observed while studying the moments in the previous section, the approach to 
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FIG. 7. The low behavior of S2cq{u,Q'^) as a function of u for various values of Q^. The 
solid curve shows = 0.5 GeV^, the short-dashed curve shows = 1 GeV^, the long-dahsed 
curve shows = 2 GeV^, and the dotted curve shows Q2 = 50 GeV^. 

scaling when using a virtual photon is slower than for the all-scalar case discussed in Q . It is 
clear that one needs to reach fairly large values of before the "resonance curve" averages 
with good accuracy to the scaling function. Indeed, with our choice of Breit-Wigner width, 
this happens only when the bumps have already disappeared, i.e. for > 5 GeV^. 

In order to illustrate this point, we have included Fig. |, where we used a value of F = 50 
MeV to smooth out the energy (5-function. The curves are more jagged than for the larger 
width, and the = 5 GeV^ curve still shows plenty of resonance structure. 

Overall, we find that the onset of local duality is definitely slower than for the all scalar 
case, which is what one expects due to the additional structure in the more realistic case 
discussed in this paper. 

IX. SUMMARY AND OUTLOOK 

We have presented a simple, quantum-mechanical model which allows us to obtain the 
qualitative features of Bloom-Gilman duality. The model assumptions we made are very 
basic: we assumed relativistic, confined valence scalar quarks, and treated the hadrons in 
the infinitely narrow resonance approximation. To simplify the situation further, we did 
not consider a three quark "nucleon" target, but a target composed of an infinitely heavy 
antiquark and a light quark. In contrast to where all particles involved in the reaction - 
electrons, photons and quarks - were considered to be scalar, we only use scalar quarks in 
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FIG. 8. The low behavior of Scq{u,Q'^) as a function of u for various values of Q^. The 
Breit-Wigner width chosen to smooth out the energy (5-functions is T = bQMeV . The solid curve 
shows = 0.5 GeV'^, the short-dashed curve shows = 1 GeV'^, the long-dashed curve shows 

= 2 GeV^, the dotted curve shows = 5 GeV'^, and the dash-dotted curve shows = 50 
GeV2. 

this paper. This makes the present model more realistic - in particular, we were able to use 
a conserved current here. However, there is still much work to be done in modelling. The 
goal of our model calculations is to gain qualitative insight into duality, its applicability and 
accuracy in various kinematic regions, not to quantitatively describe any data. In the future, 
we plan to describe more realistic situations in our model. Note that our assumptions are 
very basic and general, so that we will be able to extend our model in a straightforward 
manner. 

There are several conditions that must be fulfilled in order to see duality. In this paper, 
we put a special emphasis on three of these conditions: we demonstrated how the transition 
from coherent scattering at low to incoherent scattering at high takes place, we 
highlighted the role of relativity by considering the contributions to the Coulomb sum rule 
in a relativistic framework and a non-relativistic framework, and we gave an analytic proof 
for the equality of the scaling curve in our model and the parton model result. 

Quark-hadron duality is not only very interesting in itself, it also opens the door to very 
useful applications: duality relates the resonance region data to data from the deep inelastic 
region. If duality is understood well enough, and if the correct procedures for the averaging 
of the resonance data and the attendant errors are established, we may exploit duality 
to gather information in previously unreachable regimes. The investigation of polarized 
structure functions in the high Bjorken-x region, Xbj —>■ 1, is a major part of the experimental 
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program at Jefferson Lab Even without Icnowing details about the correct averaging 

procedures, it is clear from the experimental results and our investigation of duality that 
the conventional, sharp distinction between the "resonance region", corresponding to an 
invariant mass W < 2 GeV, and the "deep inelastic region" where W > 2 GeV, is entirely 
artificial. 

While quark-hadron duality has been investigated by theorists before, modelling duality 
is an important new step in our way to a thorough understanding of this phenomenon. In 
the literature, one often finds the phrase that duality has been explained in terms of QCD 
by DeRujula, Georgi, and Politzer What was stated in their paper is that at moderate 
Q^, the higher twist corrections to the lower moments of the structure function are small. 
The higher twist corrections arise due to initial and final state interactions of the quarks 
and gluons. Hence, the average value of the structure function at moderate is not very 
different from its value in the scaling region. While all this is true, the statement is merely 
a rephrasing of the experimentally observed fact that the resonance curve averages to the 
scaling curve in terms of the language of the operator product expansion (OPE). However, 
the operator product expansion does not explain why a certain correction is small or why 
there are cancellations - the expansion coefficients which determine this behavior are not 
predicted in the OPE. The ultimate answer to this question might come from a numerical 
solution of QCD on the lattice, but an understanding of the physical mechanism that leads to 
the small values of the expansion coefficient in the framework of a model is highly desirable. 
Also, the OPE will break down for very low Q^. Duality was experimentally observed 
to hold for as low as 0.5 GeV^ - a region where the validity of the OPE is questionable. 
In our analysis of the moments and their expansion coefficients, it became clear that a rigid 
application of the OPE at very low will inevitably lead to large, alternating expansion 
coefficients. 

The constant resonance to background ratio aspect of duality was addressed in several 
papers by Carlson and Mukhopadhyay ||23|. They used counting rules to find the de- 
pendence of the form factors of the resonances in the Breit frame, and compared them to 
the behavior of the scaling curve for large xbj and to the behavior of the background in 
the same region. From these considerations, they could explain the constant ratio, provided 
the Q^-independent coefficients of the helicity amplitudes were not anomalously small, as in 
the case of the A resonance, for which the ratio vanishes. Still, there is no explanation why 
the coefficient is small in one case and not in others, and there exist several models with 
contradictory predictions. 

The preceding observations clearly show the need for modelling. Even though one may 
obtain expansion coefficients from calculations on the lattice, an understanding of the under- 
lying physical mechanisms will most likely be gained only by considering models like ours. 
One great advantage of a purely analytical model like the one presented here is that explicit 
derivations of key quantities like the scaling function are feasible. The proof that the scaling 
function obtained for the transition from a bound quark to an excited bound quark is the 
same as the scaling function for the transition from a bound quark to a free quark was given 
here for a linear potential, which is the relativistic analog of a harmonic oscillator. It is 
desirable to extend the investigation to other types of potential, and to find a proof that 
applies to a general class of potentials. In a recent publication, IQ, numerical methods were 
applied to study the responses of a massless quark, and a disagreement between results with 
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and without final state interactions were observed. This is in contrast to our findings, and 
it is important to understand the reasons for these differences. 

The experimental data at very low still average to one single curve independent of 
their 0. However, this is not duality in the sense defined in this paper because this 
curve differs from the scaling curve. To investigate this interesting observation, one must go 
beyond the model presented here, which contains valence quarks only, and therefore must 
produce a valence-like shape. However, introducing sea quarks and modeling the decay of 
the excited resonances, along with the corresponding non-resonant production mechanisms 
from sea quark pairs, might shed considerable light on this issue. 
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APPENDIX A: DUALITY IN THE EXCITATION FORM FACTORS 

In this appendix, we proceed to study duality in the excitation form factor FQ^^{(f). 
While this duality is not directly related to an observable like the structure functions or 
response functions, it exhibits duality very clearly. Duality in the form factors has recently 
received some attention in 

The duality prediction is that it should not matter if we describe the process in question 
in a perturbative QCD picture, involving only a free quark, or in a hadronic picture with 
resonances. The form factor for a hypothetical free quark is just 1, as it does not have 
any structure. In the hadronic picture, we have inclusive electron scattering where we can 
excite all resonances - as the final hadronic state is not observed, we have to sum over the 
resonances incoherently. So we have to compare l-^o,Af(^^)P to 1. Fig. ^ shows single 
form factors for the lowest resonances, the elastic peak and inelastic excitations up to N = 
5. All form factors look qualitatively similar, except for the N=0 elastic form factor, which 
starts at 1 for |g| =0. In general, the form factors increase in width and decrease in height 
with increasing N. 

Using our previous expression for the form factor, eq. (|13D, we find for the sum up to a 
certain value Nmnr 
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FIG. 9. The excitation form factor squared, |i*bAf(9*)|^ for the lowest excitations, N = to 
N = 5, and /? = 0.4 GeV. 

as mentioned when discussing the Coulomb sum rule. However, we are limited in the maxi- 
mal value of Nmax not by technical problems, but by a physical constraint: we are considering 
electron scattering, i.e. space-like kinematics, and therefore we must fulfill the condition 

>0^\q\> iy = En- Eo, (A3) 

so that for fixed three-momentum transfer we find a limit on the value of N. The form 



factor sums are shown in Fig. |10| for various values of the oscillator parameter /3. Larger 
values of (3 indicate a stronger binding. 

The spiky character of the curves stems from the fact that the form factors for the N-th 
state are only allowed to contribute to the sum if |g| > un, where '■= E^ — Eq. Therefore, 
the sum jumps up whenever another threshold is crossed. This effect can be observed best 
for the strongest binding, j3 = 0.6 GeV, as the gaps between the energy levels are largest in 
this case. With increasing value of the three-momentum transfer, more and more resonance 
states can be excited and contribute to the sum: the spikes subside and the average value 
of the sum gets fairly close to 1. The curve for the weakest binding, (3 = 0.2 GeV, becomes 
almost smooth and takes on a value of 0.9994, i.e. duality is violated by less than 0.06 %. 
For jS = 0.4 GeV, the sum reaches 0.995, so duality is violated only by 0.5 %. Even for the 
strongest binding, /3 = 0.6 GeV, the duality prediction is fulfilled within 1%. Here, we see a 
typical feature of duality, the need for many resonances to contribute in order to reproduce 
the behavior of a free quark. At low where only a few resonances can contribute, the 
deviations from 1 are larger. The fact that duality is fulfilled best for weak binding is what 
we expect: a quark that is bound very lightly and then receives a hard momentum transfer 
behaves essentially as if it was free. If the binding gets stronger, the situation gets more 
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FIG. 10. The sum over the excitation form factors squared, for /? = 0.2 GeV (soHd hue), 
P = 0.4 GeV (dashed Une), and /? = 0.6 GeV (dash-dotted hue). 

non-perturbative, and duality does not work as well. 

The duality as seen in the form factors is reminiscent of duality in the decay rates of the 



semileptonic decay of heavy quarks [11,27|. There, in the limit of infinite masses of the b 
and c quarks, the loss of strength in the elastic channel is compensated for by the increase 
in the inelastic decay channels. Once one considers heavy, but not infinitely heavy, quark 
masses, one obtains a jagged structure, with peaks getting close to the free quark limit, quite 
similar to what we observe when considering the Coulomb Sum Rule and the excitation form 
factors. 

Let us consider the number of resonances needed in more detail, so that we can draw 
further conclusions on the kind of duality we are observing here. In the calculations presented 



in Fig. |T0|, we summed up to the highest allowed N, Nmax, which is quite large in general. 
In Fig. |ll], we present the full curve, and three curves where we summed over a limited 
number of resonances only, namely from N = 10 to N = 40 (dash-dotted line), N = 10 to 
N = 30 (dotted line) and N = 20 to N = 40 (dashed line). One can see that for a small 
interval in three-momentum transfer of g = 2.6 GeV - 3.0 GeV, it is sufficient to include 
only resonances from N = 10 to N = 40 in order to reproduce the full, unrestricted curve. 
One also sees from the dotted line that for lower three-momentum transfer, the inclusion 
of just the twenty resonances from N = 10 to N = 30 suffices to get close to the full curve 
q = 2.6 GeV, while the same number of resonances does not suffice to approximate the full 
curve at a slightly higher value of |g| (see the dashed line). 

So, while it is clear that the "degrees of freedom duality" holds very nicely over the whole 
kinematic range, we see that duality - the truncated version - does not hold as well: we do 
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FIG. 11. The sum over the excitation form factors squared |-Fo_Ar(g^)p for /? = 0.4 GeV. The 
sohd Hue shows the sum over all allowed N, the dash-dotted line shows the sum from N = 10 to N 
= 40, the dotted line shows the sum from N = 10 to N = 30, and the dashed line shows the sum 
from N = 20 to N = 40. 

need a certain number of resonances to obtain duality in a limited kinematic interval, and 
this number increases when we increase the three-momentum transfer. 

APPENDIX B: THE ROLE OF RELATIVITY 

We have stressed the importance of relativity before, and while it is quite obvious that 
one needs to include it for GeV momentum transfers to light quarks, it is instructive to see 
how relativity works for the form factors. In order to illustrate this point, we show the sum 
of the excitation form factors squared calculated for the non-relativistic harmonic oscillator 
potential in Fig. 0. As mentioned in section |l[ the wave functions, and therefore the form 
factors, are the same. The difference lies in the energy spectrum. 

Obviously, duality in the non-relativistic case does not work at all: the curves start out 
at 1, as the elastic form factor for |g| = is 1, but then fall off immediately. Whenever a 
new threshold opens, the additional contribution is not sufficient to compensate the fall-off 
of the other form factors: they do not contribute at their maximum value, but only with 
the "high side of the peak, where the form factor drops quickly, see Fig. |^. In the 
non-relativistic case, the spacing between the energy levels is wider than in the relativistic 
case, where the levels shrink together. This means that considerably fewer resonances are 
allowed to contribute at the same, fixed and therefore the resulting sum is much smaller 
and duality is violated. To clarify this point, we show a comparison of energy levels in 
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FIG. 12. The sum over the excitation form factors squared for the non-relativistic case, for 
P = 0.2 GeV (soHd hne), P = 0.4 GeV (dashed Hne), and /? = 0.6 GeV (dash-dotted hne). Note 
that the scale differs from the scale in Fig. |l^ 



Fig. |T3. 

In the left panel, we show the energy transfers u for the first eleven energy levels (elastic 
and the first ten resonances) of a relativistic harmonic oscillator with f3 = 0.4 GeV. The 
diagonal line marks the photopoint, \q\ = u. This means that for a given all the energy 
transfers below that line are in the space-like region and therefore allowed. E.g., in the 
relativistic case for |g| = 1 GeV, nine resonances (counting the elastic) can be excited. The 
right panel of Fig. |1^ shows the energy transfers u for the energy levels of a non-relativistic 
harmonic oscillator with (3 = 0.25 GeV. We chose the (3 for the non-relativistic oscillator in 
order to reproduce the energy splitting between the ground state and the first excited state 
of the relativistic case. One clearly sees that fewer resonances can contribute here, e.g. only 
six compared to nine in the relativistic case at |g| = 1 GeV. The discrepancy grows larger 
for higher as the relativistic energy levels move closer together, while the non-relativistic 
ones are equally spaced. 

In conclusion, we have seen that the relativistic description is necessary to ensure a 
correct treatment of the phase space. Only with a proper relativistic phase space do we see 
duality in the excitation form factors. This was already clear from our discussion or the 
Coulomb Sum Rule in the main body of the paper. 

Mathematically, degrees of freedom duality in the excitation form factor means that 
one can expand a plane wave (free quark) in a set of Hermite polynomials (bound quark), 
provided one uses a sufficiently large number of basis states. Any other set of ortho-normal 
polynomials would also work. 
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FIG. 13. The first eleven energy transfers v for the relativistic case with /3 = 0.4 GeV (left 
panel) and for the non-relativistic case with /3 = 0.25 GeV (right panel). The diagonal line in both 
panels indicates the photopoint, i.e. |g| = v. All energy transfers below this line are allowed in the 
space-like region. 
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